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Abstract 

We present the calculations of the complete next-to-leading order (NLO) inclusive total cross 
sections for the associated production processes pp UXk + X in the Minimal Super symmetric 
Standard Model at the CERN Large Hadron Collider. Our calculations show that the total cross 
sections for the tiXi pi'oduction for the lighter top squark masses in the region 100 GeV < < 
160 GeV can reach 1 pb in the favorable parameter space allowed by the current precise experiments, 
and in other cases the total cross sections generally vary from 10 fb to several hundred fb except 
both > 500 GeV and the ^2X2 production channel. Moreover, we find that the NLO QCD 
corrections in general enhance the leading order total cross sections significantly, and vastly reduce 
the dependence of the total cross sections on the renormalization/factorization scale, which leads 
to increased confidence in predictions based on these results. 
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I. INTRODUCTION 



The CERN Large Hadron Collider (LHC), with = 14 TeV and a luminosity of 100 
fb~^ per year |l^, will be a wonderful machine for discovering new physics. In so many new 
physical models, the Minimal Supersymmetric Standard Model (MSSM) P| is one of the 
most attractive models for the theorists and the high energy experimenters, and searching 
for supersymmetric (SUSY) particles, as a direct experimental evidence, is one of the prime 
objectives at the LHC. Therefore, a good understanding of theoretical predictions of the 
cross sections for the production of the SUSY particles is important. 

The cross sections for the production of squarks and gluinos were calculated at the Born 
level already ^any years ago ^j. To date, tie productions of gluinos and squarks jj, 0], 
top squarks sleptons US] and gauginos 3| at the hadron colliders have been studied. 
Besides the leading order (LO) results, these calculations also have included the next-to- 
leading order (NLO) corrections to improve the accuracy of the theoretical predictions and 
reduce the uncertainty from the dependence of results on renormalization/factorization scale. 
It was recently pointed out in Ref. that the associated production of a gaugino (x) with a 
gluino {g) or with a squark (g) is potentially a very important production mechanism, since 
the mass spectrum favors much lighter masses for the color-neutral charginos and neutralinos 
than for the colored squarks and gluinos in popular models of SUSY breaking. And a detailed 
NLO QCD calculation of the associated production of a gluino and a gaugino at the Tevatron 
and the LHC has been given in Ref. j^. However, the processes of the associated gaugino 
and squark productions have not been calculated at LO and NLO yet. 

In this paper, we present the complete NLO QCD (including SUSY QCD) calculation for 
the cross sections of the associated production of top squarks (stops) and charginos at the 
LHC. Similar to pp ^ gb ^ tH^ 10], which is expected to be a dominate process for the 
charged Higgs boson production at the LHC, the associated production pp gh iiXk 
may be also the dominate process for single top squark or chargino production at the LHC. 
This is due to the following reasons: first, the large top quark mass in stop mass matrix can 
lead to strong mixing, and induce large mass difference between the lighter mass eigenstate 
and the heavier one, which means that the phase space for the lighter stop will be great 
and benefit its production; second, besides containing a strong QCD coupling between the 
incoming partons, this process also includes an enhanced effects from the Yukawa coupling 
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in the vertex b — ti — Xk final states. For simplicity, in this paper we neglect the bottom 

quark mass except in the Yukawa coupling. Such approximations are valid in all diagrams, 
in which the bottom quark appears as an initial state parton, according to the simplified 



Aivazis-Collins-Olness-Tung (ACOT) scheme However, it was pointed out in Ref 
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that the approximations of the hard process kinematics and the introduction of conventional 
bottom quark densities will give rise to sizable bottom quark mass and kinematical phase 
space effects, and may overestimate the inclusive cross section, for exarnple, in the processes 
of 66 — >■ if and bg —>■ bH [1^. Very recently, it is shown in Ref. that the bottom 
parton approach is still valid if we choose the factorization scale below the average final 
state mass: fif ~ Cmav = C{mi. + m--)/2 with C ~ (1/4, 1/3). Thus, in this paper, we 
choose /i/ = ?Tiav/3 when we use the bottom parton approximation. 

The arrangement of this paper is as follows. In Sec. II we show the LO results and define 
the notations. In Sec. HI we present the details of the calculations of both the virtual and 
real parts of the NLO QCD corrections. In Sec. IV by a detailed numerical analysis we 
present the predictions for inclusive and differential cross sections at the LHC The lengthy 
analytic expressions are summarized in Appendixes A and B. 

II. LEADING ORDER PRODUCTION OF STOPS AND CHARGINOS 

The partonic process of the LO associated production of stops and charginos is 
g{pa)b{ph) UipilXkiPi), and the related Feynman diagrams are shown in Fig. [H In 
order to simplify our expressions, we introduce the following Mandelstam variables: 

s = {Pa+ Pbf, sa = s- ml - ml- , 
t={Pa-pif, ti=t-ml, t2=t-m\-, 

u={Pa-P2f, Ui = u-ml, U2=u-m\-. (2.1) 

After the n-dimensional phase space integration, the LO partonic differential cross sections 
are given by (n = 4 — 2e) 

J2-B q t2U2-Sm^,- 

- ' ' '''-rQ{t2U2-sml-)Q[s-{mi^ + m--f] 



dt2du2 s2r(l - e) ' /i2s ' ' ^fe ' ' ' '^fc 



x6{s + t + u-ml- m'^.)\Mi\ (2.2) 
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with 



12 l-e 



(I4r + i^-.i')( 



t2 + {l-e)u2 ss^-ui{u2 + 2t2) 2t2m| 



sU 



+ 



(2.3) 



where = (47r) and the indices {i,k) label the outgoing particles {ti,Xk)- The scale 
parameter fir is introduced to provide the correct mass dimension for the coupling constant 
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in n-dimensions. \M^\ is the LO squared matrix element, which has been summed the 
colors and spins of the outgoing particles, and averaged over the colors and spins of the 
incoming ones. Ijj^ and fc*^ are the left- and right- handed coupling constants of the vertex 
h — U — Xk 1 respectively, and are defined together with and k\f, in the vertex t — hi — xt-i 
which is involved in the virtual corrections, as follows: 



with 



nit 



+ 



YbRi2Uk2 
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ut 



Ub 



nib 



YbR\iUk2 
YtRiVk2 



h 
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(2.4) 



(2.5) 



V2miv sin /3 ' V2niiY cos (3 

Here the angle f3 is defined by tan/? = V2/V1, the ratio of the vacuum expectation values 
of the two Higgs doublets. Matrices U and V are the chargino transformation matrices 
from interaction to mass eigenstates defined in Ref. 16|. Matrix R^ [q = t,b) is defined to 
transform the current eigenstates and qR to the mass eigenstates qi and q2- 



R" 




R'i 



cos 6q sin 6q 
— sin 9q cos 9q 



(2.6) 



with < < vr by convention. Correspondingly, the mass eigenvalues mq^ and niq^ (with 



TT^qi < ""^52) given by 



/ 2 



m4 



Rm?{R'^ 



/ 2 

aqniq 



yaqniq m| 



(2.7) 



with 



M| + m| cos 2/5(J|^ - Cg sin^ Ow) 



Cg = Ag — /x{cot/5, tan/3} 



(2.8) 
(2.9) 
(2.10) 
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for {up, down} type squarks. Here M| is the squark mass matrix. Mq qj^ and At^b are soft 
SUSY-breaking parameters and yU is the Higgs mixing parameter in the superpotential. 
and Cq are the third component of the weak isospin and the electric charge of the quark g, 
respectively. 

The LO total cross section at the LHC is obtained by convoluting the partonic cross 
section with the parton distribution functions (PDFs) Gb,g/p in the proton: 

'^fk = J dxidx2[Gb/p{xi, fif)Gg/p{x2, fXf) + (xi ^ X2)]ajl, (2.11) 

where /i/ is the factorization scale. 



III. NEXT-TO-LEADING ORDER CALCULATIONS 

The NLO contributions to the associated production of stops and charginos can be sep- 
arated into the virtual corrections arising from loop diagrams of colored particles and the 
real corrections arising from the radiations of a real gluon or a massless (anti) quark. 



A. Virtual Corrections 

The virtual corrections to gb tiXk arise from the Feynman diagrams shown in Fig. |21 
and 01 which consist of self-energy, vertex and box diagrams. We carried out the calcu- 
lation in t'Hooft-Feynman gauge and used the dimensional regularization in = 4 — 2e 
dimensions to regularize the ultraviolet (UV), soft infrared and collinear divergences in the 
virtual loop corrections. However, this method violates the supersymmetry. In order to 
restore the supersymmetry the simplest procedure is through finite shifts in the quark- 



squark-gluino(chargino) couplings [17[: 

n 4 

e = e(l-^Cp), 
Svr 

hh = Kb{l-^-^CF) (3.1) 

with = 3 and Cp = 4/3. Since we have not such a coupling as ^ — g — g at the tree-level, 
the first shift is not used in our calculations, and we just consider the latter two shifts. As 
for the Dirac matrix 75, we deal with it using the "naive" scheme, in which the 75-matrix 



anticommutes with the other 7^,-matrices. This is a legitimate procedure at the one-loop 



level for anomaly-free theories 



The virtual corrections to the differential cross section can be expressed as 



dt2du2 s^T{l - e) ' fi^^s ' ' ' ' 

x5(s + t + u - m| - m|-) 2 Re(M^M|*), (3.2) 

where M^^ is the renormalized amplitude, which can be separated into two parts: 

= MT""" + (3.3) 

Here Mj'^"^*^" contains the self-energy, vertex and box corrections, and Mj^^" is the corre- 
sponding counterterm. Moreover, M™*"^" can be written as 

Munren = j:M:,+ ±M^^<'\ (3.4) 
a=a /3=a 

where a and (3 denote the corresponding diagram indexes in Fig. |2l and Fig. IHl respectively. 
We express each term further as 

Mt, = Csf:f^Mu (3.5) 
1=1 

m;^°"(^) = c,f;/^"-(^)M;, (3.6) 

1=1 

where Cg = igg^ /{IGTr"^), ff and fj'"^^'^^ are the form factors, which are given explicitly in 
Appendix A, and Mi are the standard matrix elements defined as 

Mi,2 = V^{P2)C~'T^ APa)PL,R u{p,), 
Ms,4 = V^{p2)C-^T^ A MPa)PL,R u{pb), 

M5,6 = v''{p2)C-^T-Pl,r u{p,)p, ■ e{pa), 

M7,8 = V^{P2)C-^T-Pl,R U{p,)p2 ■ 6{pa), 
M9,io = V^{p2)C-'T^ ^aPL,R u{pb)pb " e{pa), 

Mn,12 = V^ip2)C-'T'^ ^aPL,R w(Pb)P2 " ^(Pa). (3.7) 

Here C is the charge conjugation operator, and T" is the SU{3) color matrix. 
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Mj™" is separated into Mj™"*-*-* and Mj™"*-*-*, i.e. the counterterms for s and t channels, 
respectively, which are given by (j 7^ i) 



j^conis) ^ _!^{[4a + li^sz^, + 54](M3 + 2M5) + [4a + + 54](M4 + 2M,)}, 

+ "-^ + ^^^]{Me-Ms)}, 



ti ti ti t — mf ti ti 



with 



Silk = Rl2VkW + >^*K2(^*2^ + Rli^o' 



nit - • ' 



54 = n^2(4^-4'5^')- 

where Smf , Srrig (g = t,b), SZf,, SZg, SZ^-, SZjj and S9f are the renormalization constants, 
which are fixed by the on-shell renormalization scheme [l9| , and can be written as (assuming 

S^i = -^C^{4m? 5o(m|, m| , 0) - Ao(m? ) + 4Ao(m?) + 4[m? 5i + m^i^o 
-2mtmgR\^R\^BQ]{m\,m], m^)}, 

^ = -^Cp{(45o + 250K,mJ,0) - 1 + Y.[Bi - ^ sin2^,-(-l)^i?o]K, m|, m? )}, 

SZ,= ^CpB,{0,mlml), 
a 2 

SZg = -^[-{Bo + 2m]B',) (0, m?, m?) + 2(So + 2m2s^) (0, m^, m?) 
+^E(^^o-m?fi;)(0,m?,m?)-^], 

= ^CH(5o + 2m? 5;)(m|, m|, 0) - [5o + K " " "^?)^o 



27r 

+4.mtm~gR\^R\^B'Q] {ml ,ml,ml)}, 



- ReSUm? 
5Z' ' ' 



^-^ mf — mf ' 



Re[Sl,(m|) + Sl,(m|)] 
2(m|-m|) 
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where Bq = dBo/dp'^, Aq and Bi are the Passarino-Veltman one-point and two-point in- 
tegrals, respectively, which are defined similar to Ref. Q except that we take internal 
masses squared as arguments. Here the counterterm 66i for the stop mixing angle is fixed 



as Ref. 



2l|, and 



a 2 

= ^-{sin4%[Ao(m|) - Ao(m|)] + Sm^m^ cos 2^,-So(m? , m|, m^)}. 



In SZ^-, i?Q(m?, m?, 0) contains infrared divergences, and can be expressed as 



B'.iml m|, 0) = -^(^)T(1 + + 2). 
Moreover, the QCD coupling constant Qs is renormalized in the modified minimal subtraction 



(MS) scheme except that the divergences associated with the top quark and colored SUSY 
particle loops are subtracted at zero momentum l^l. Assuming that the scalar partners of 
the Hf = 5 light quark fiavors have a common mass rrig, we obtain 

2\ /3 1 AT rm"^ ^ _ .m2 1 m? 1 m? 1 »t^2 



^9s _ /3ol , N m\ n.m\ 1 1 1 



+ — ln^ + ^ln^ + — ln^ + — ln^ + -ln ^1,(3.8) 

which agrees with one given in Ref. jj]. Here 1/e = l/ejyy — 7b -|- ln(47r), and 

/3o = (yiV - \n^) + [-^(iV + 1) - + 1)] ^ /^^^ + /3o^, (3.9) 

where (5^ includes the contributions from the gluon and the quarks except the top quark, 
and contains the contributions from the top quark and all colored SUSY particles. The 
evolution of is determined by (5^. After renormalization, M^^ is UV-finite, but it still 
contains the infrared (IR) divergences: 



M^lzH = f^r(l + + (3.10) 

zTT see 

where 



A\ = -f, (3.11) 



aV 43 4 3 -h 1 -ui 

A[ = — - - -In— + 31n— - -In^-. (3.12) 
o 3 rrif rrif 3 rrif 

LA Ln Li 



Here the infrared divergences include the soft infrared divergences and the collinear infrared 
divergences. The soft infrared divergences can be cancelled by adding the real corrections, 
and the remaining collinear infrared divergences can be absorbed into the redefinition of 
PDF 0, 

which will be discussed in the following subsections. 
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B. Real Gluon Emission 



The Feynman diagrams of the real gluon emission process g{pa)b{pb) — > ii{pi)xl (P2) + 
qIps) are shown in Fig.4. Using the notations analogous to ones defined in Ref. |2J], we 
express our results in terms of the following invariants: 

S = {Pa+ Pbf, S5 = (Pl + J^2)^ S5A = S5 - A, 



S3 = {P2 + Psf - rn\, S4 = {pi + Pzf - m|, Ssa = S3 - Af, 
t={Pb-P2f, t' = iPb-P3f, ti = t-ml, t2 = t-ml-, 

U=iPa-P2f, U=iPa-P3f, Ui = U-m]_, U2 = U - IV? 

U6 = {Pb - Pif - m\, Ut = {pa - pif - m\, (3.13) 

where A^ = mf + m?_ and Aj = — m? . 

The phase space integration for the real gluon emission will produce infrared singularities, 
which can be either soft or collinear and can be conveniently isolated by slicing the phase 
space into different regions defined by suitable cutoffs. In this paper, we use the two cutoff 
phase space slicing method j3] to decompose the three-body phase space into three regions. 

First, the phase space can be separated into two regions by an arbitrary small soft cutoff 
5s, according to whether the energy of the emitted gluon is soft, i.e. < 5sy^/2, or hard, 
i.e. > 6s\^/2. Correspondingly, the partonic real cross section can be written as 

^fk = ^i + ^fk, (3.14) 

where and afl are the contributions from the soft and hard regions, respectively, af^ 
contains all the soft infrared divergences, which can be explicitly obtained after the integra- 
tion over the phase space of the emitted gluon. Second, to isolate the remaining collinear 
divergences from afl, we should introduce another arbitrary small cutoff, called collinear 
cutoff 6c, to further split the hard gluon phase space into two regions, according to whether 
the Mandelstam variables satisfy the collinear condition < —t' (— m') < ScS or not. Then 
we have 

^.^ = '^.r + ^F, (3.15) 

where the hard collinear part a^'" contains the collinear divergences, which can be explicitly 
obtained after the integration over the phase space of the emitted gluon. And the hard non- 
coUinear part is finite and can be numerically computed using standard Monte-Carlo 



integration techniques 26|, and can be written as 



rf^r = ^|Mif'|V3. (3.16) 

Here dTr, is the hard non-coUinear region of the three-body phase space, and the exphcit 
„^ of Ml ca„ be fou„d . AppeM. B. I„ o J. I avo.d .n.od.c.. eln. 
ghost hnes while summing over the gluon hehcities, we hmit ourselves to the sum over the 
physical polarizations of the gluons j2J], i.e. 

where the index i (=1,2) labels the two external gluons, and 7^ ki is an arbitrary vector. 
This polarization sum obeys the transversality relations 

k^.P^" = P'^^h, = n,,P^'' = P'^^n,, = 0. (3.18) 

In the next two subsections, we will discuss in detail the soft and hard coUinear gluon 
emission. 



1. Soft gluon emission 

In the limit that the energy of the emitted gluon becomes small, i.e. < Ssy/s/2, 

-2 



the matrix element squared |Mj^| can be simply factorized into the Born matrix element 
squared times an eikonal factor ^eik- 

\Mg{gh^Ux,+9)f iAnasfi'/)jMi\\,,,, (3.19) 

where the eikonal factor $ejfc is given by 

= N^^!^l^^— + N ^" " - C, ""^^ - ^ .(3.20) 

{Pa ■ P3)iPl ■ P3) {Pb ■ PzliPa ■ P-i) {Pl-PzY N {pi ■ P3){pb ■ Ps)' 

Moreover, the phase space in the soft limit can also be factorized as 

dr3{gb UXk + 9) dT2{gb Ux^)dS, (3.21) 
where dS is the integration over the phase space of the soft gluon, which is given by 

2(27r) 



= ^3 2e / rf^3^3"''^^2-2.. (3.22) 
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Then the parton level cross section in the soft region can be expressed as 

af, = {Airasfil') J rfFsp^' J dS^e^k. (3.23) 



Using the approach of Ref. 3|, 
Eq. ()3.23|) becomes 



after the integration over the soft gluon phase space, 



*--*"l2^rra(^)'K7^ + f ^-^S) (3.24) 



with 

13 



A 



2 



3 



A? = -2A^ln5, + ^ + ^ln^-31n^ + ^ln^, 
6 6 rrij rrif 3 

ti t'l 

= 2A^ In^ - 2AI In (5, + {Cp^ - ]) In + ^ In^ ^i^^jl + Li2[l + 



ti 6 ui 2 7-/3 s(/3-7) 

-2Li.i ''+y-^% . 



where 7 = (s + m| — "m^-)/ (2s) and /? = ^72 _ mf /s. 



2. Hard coUinear gluon emission 

In the hard collinear region, i.e. E3 > 6sy/s/2 and < —t' {—u') < 6cS, the emitted 
hard gluon is collinear to one of the incoming partons. As a consequence of the factorization 



theorems [27|, the squared matrix element for gb —* tiX^ + 9 can be factorized into the 
product of the Bor n sq uared matrix element and the Altarelli-Parisi splitting function for 



hg and g ^ gg 



3, i.e. 



\Mg{gh - Ux, + 9)\" (4vra./.-)|M^r(^^^ + ^^^), (3.25) 

where z denotes the fraction of incoming parton b{g)'s momentum carried by parton b{g) 
with the emitted gluon taking a fraction {1 — z), and Pij{z, e) are the unregulated splitting 
functions in n = 4 — 2e dimensions for < z < 1, which can be related to the usual 
Altarelli-Parisi splitting kernels j3| as Pij{z,e) = Pij{z) + ePlj{z), explicitly 

P,,iz) = 2N[-^ + 1^ + ^(1 _ ^)], P'jz) = 0. (3.26) 
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Moreover, the three-body phase space also can be factorized in the coUinear hmit, and, for 



example, in the limit < —t' < 6cS it has the following form j25|: 



Note that the two-body phase space should be evaluated at a squared parton-parton energy 
of zs. Then the three-body cross section in the hard collinear region is given by 2^ 

dz 1 — z 

+Pgg{z,e)Gg/p{xi/z)Gb/p{x2) + (xi ^ X2)] — {^^)~^dxidx2, (3.28) 
where Gb{g)/p{x) is the bare PDF. 

C. Final States with an Additional Massless Quark 

In addition to the real gluon emission, other real emission corrections to the inclusive cross 
section for pp iiXk NLO involve the processes with an additional massless (anti)quark 
in the final states, as shown in Fig.5 {q = u, d, s, c): 

9{Pa) + g{Pb) ^ ti(pi) + Xk{P2) + KP3), (3.29) 

?/^(Pa) + h{pb) u{pi) + Xk{P2) + q/qips), (3.30) 

b/bipb) + bipa) ^ iM + xM + VKps), (3.31) 
q{Pa) + q{Pb) Upi) + Xk{P2) + Kpz)- (3.32) 

Since the contributions from the processes (j3.29|) - (j3.3ip contain the initial state collinear 
singularities, we also need to use the two cutoff phase space slicing method [25]. However, 
we only split the phase space into two regions, because there are no soft divergences here. 



Thus, according to the approach shown in Ref. 



the cross sections for the processes with 



an additional massless (anti)quark in the final states can be expressed as (g = u, d, s, c, b) 
dc^ik'^ = ^ik{oi(^ iiXk + X)[G a/p{xi)G i3/p{x2) + (a;i ^ X2)]dxidx2 

(",/3) 

+d^f^i^ l[lZ2e) (^)^K-7)^c-ln,(^, e)Gg,p{x,/z)Ggip{x2) 

+ Pga{z,e)Go,/p{xi/ z)Gh/p{x2) + {xi ^ X2)] — ( )~^dxidx2, (3.33) 
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where 



1 _L fl _ zV 



z 

P,,iz) = ^[z' + (l- zf], P;(z) = -zil - z). (3.34) 

The first term in Eq. ()3.33|l represents the non-collinear cross sections for the four processes, 
and hke Eq. ()3.1(i|l each can be written in the form: 

1 



rfag = -|Mf/| dla, (3.35) 

where a and /3 denote the incoming partons in the partonic processes, and dV^ is the three 
body phase space in the non-colhnear region. The explicit expressions of |Mj^ | can be 
found in Appendix B. The second term in Eq. ()3.33p represents the colhnear singular cross 
sections. 

^ ^ ^ — 

Moreover, for the subprocesses gg/qq — > tit- — > tiX^b {q = u,d, s,c,b), assuming m^. > 
m^-, the stop momentum can approach the stop mass shell, which will lead to singularity 
arising from the stop propagator. Following the analysis shown in Ref. |^, this problem 
can easily be solved by introducing the non-zero stop widths F^. and regularizing in this 
way the higher-order amplitudes. However, these on-shell stop contributions are already 
accounted for by the LO stop pair production with a subsequent decay into a chargino and 
a bottom quark, and thus should not be considered as a genuine higher order correction to 
the associated production of top squarks and charginos. Therefore, to avoid double counting, 
these pole contributions will be subtracted in our numerical calculations below in the same 



.a. 



way as shown in Appendix B of Ref. 



D. Mass factorization 

As mentioned above, after adding the renormalized virtual corrections and the real cor- 
rections, the partonic cross sections still contain the collinear divergences, which can be 
absorbed into the redefinition of the PDF at NLO, in general called mass factorization [2^. 
This procedure in practice means that first we convolute the partonic cross section with the 
bare PDF Ga/p{x), and then use the renormalized PDF Ga/p{x, fif) to replace Ga/p{x). In 
the MS convention, the scale dependent PDF Ga/p{x, ^f) is given by [2^ 

G^/p{x,^^f) = Ga/p{x)+Y.i--J[^ ^^^I^l^ i^r] [ ^Pa,{z)G,/p{x/z). (3.36) 
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This replacement will produce a collinear singular counterterm, which is combined with the 
hard collinear contributions to result in, as the definition in Ref. M, the 0{a.) expression 
for the remaining collinear contribution: 

a=b,g ^ 

+ ^ X2)}dxidx2, (3.37) 



where 



AT{q ^ qg) = CFi2\n5s + 3/2), (3.38) 
ATig ^ gg) = 2N\n6s + {IIN - 2n/)/6, (3.39) 



s 



A- = Arin(-), (3.40) 
Ga/p{x,fif) =Y, —Gf3/p{x/y,fif)Pa/3{y) (3.41) 



with 



PaM = Pap ln(4^4) - KM- (3-42) 

y fj'f 



Finally, the NLO total cross section for pp tiXk factorization scheme is 

(^ik^° = j {dxidx2[Gb/p{xi, fif)Gg/p{x2, /i/) + xi ^ X2]{afk + a(k + af^ + ^,^^) + dalf] 
+ dxidx2[Gc,/p{xi,fif)Gf3/p{x2,fif) + {xi'^X2)]af,,{a(3-*iiXk+X). (3.43) 

Note that the above expression contains no singularities since + ^2 = and + Al + 
Arib ^ bg) + Arig gg) = 0. 

E. Differential Cross section in the Transverse Momentum 

In this subsection we present the differential cross section in the transverse momentum 
Pt of the charginos. Using the notations defined in Ref. [J], the differential distribution with 
respect to pr and y for the processes 



p{Pa) +p{Pb) ^ k{pi) + Xk{P2)[+g{P3)/q{P3)/q{P3)] (3.44) 
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is given by 



2ptSY^ I _dxi f_dx2XiGa/p{xi,fif)x2Gi3/p{x2,fif) ,,^"^ , (3.45) 

apTdy Q, p Jx^ Jx^ CLt2CLU2 

where is the total center-of-mass energy of the colhder, and 

2 T2U2 2 1 , / ^2 N 

—T2 — m-_ + mf —X1U2 — mi^ + mf 
= ^2 = ^ (3.46) 

O + t/2 Xib + i2 

with T2 = (Pf, — ^2)^ Tni- and f/2 = (Pa ^ P2Y ~ fn"^--- The hmits of integral over y and 
Pt are 

- 2/"'"'=(pt) < 2/ < 2/""'^'(Pt), < pt < Pt^^, (3.47) 



with 



S" + m,- — 

^-'^-(py) = arccosh( 



k 



1 



^rnax ^ (S + m^._ - m? )2 - Am^._S. 



(3.48) 



IV. NUMERICAL RESULTS AND CONCLUSION 



We now present the numerical results for total and differential cross sections for the 
associated production of top squarks and charginos at the LHC In our numerical calculations 
the Standard Model (SM) parameters were taken to be aewif^w) = 1/128, m^y = 80.419 
GeV, mz = 91.1882 GeV, and nit = 174.3 GeV js^. We use the two-loop evolution of as{Q) 
MMz) = 0.118) and CTEQ6M (CTEQ6L) PDFs throughout the calculations of 
the NLO (LO) cross sections. Moreover, in order to improve the perturbative calculations, 
we take the running mass mi,{Q) evaluated by the NLO formula [3^: 



mb{Q) = UG{Q,mt)U5{mt,mb)mb{mh) (4.1) 
with mb{mb) = 4.25 GeV [s^. The evolution factor Uf is 

as{Qi) 47r 

dif) = -Jl- JU) = 8982 - 504/ + 40/^ 

33-2/' 3(33 -2/)2 ' ^ ' ' 

15 



In addition, in order to improve the perturbation calculations, especially for large tan P, we 
make the following replacement in the tree- level couplings j33|: 



2a, , , , , , h 



2 



= —^Mgij,tan(3I{mi^,mi^,Mg) + — ^/i^ tan/3/(mj^, mf^, /i) 



^/iM2tan/3^[(4)2/K^,M2,/i) + ^(4)2JK,M2,/x)] (4.4) 

2=1 



with 

b, c) = (,._^.)(^/,.)(,._,.) («^^^ log ^ + h-c- log ^ + c-a- log ^). (4.5) 

And, it is necessary for avoiding double counting to subtract these (SUSY-)QCD corrections 
from the renormalization constant 5mi, in the following numerical calculations. The MSSM 
parameters are determined as follows: 

(i) For the parameters Mi, M2 and /i in the chargino and neutralino matrices, we take 
M2 and /i as the input parameters, and assuming gaugino mass unification we take Mi = 
(5/3) tan2 ^i^Ma and M~g = (a,(Mg)/a2)M2 H. 

(ii) For the parameters in squark mass matrices, we assume Mq = Mjy = Mj^ and 
At = Ah = 300 GeV to simplify the calculations. Actually, the numerical results are not 
very sensitive to At(^i,y 

In our calculations, except in Fig. 10, we always choose the renormalization scale fir = m&v, 
and the factorization scale /i/ is fixed to mav/3 instead of mentioned in Sec. I. 

Before we discuss in detail the numerical results for associated production of top squarks 
and charginos at the LHC, in Fig.6 we first show that the following NLO QCD predictions 
do not depend on the arbitrary cutoffs 5s and 5c in the two cutoff phase space slicing method. 
Here we take = -200 GeV, M2 = 300 GeV, m^-^ = 250 GeV, tan/3 = 30, and 5^ = 5,/100. 
(Tother includes the contributions from the Born cross section and the virtual corrections, 
which are cutoff-independent. We can see that the soft plus hard collinear contributions 
and the hard non-coUinear contributions depend strongly on the cutoffs, and especially for 
the small cutoffs {5s < 10"'^) their magnitudes can be ten times larger than the LO total 
cross section. However, the two contributions {<Jsoft + o'hard/cou and (Thard/non-coii) balance 
each other very well, and the final result ct^lo is independent of the cutoffs and keeps 0.23 
pb. Therefore, we take 5s = 10~^ and 5c = 5^/100 in the numerical calculations of Fig. 7-13. 
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Fig. 7 shows the dependence of the LO and NLO predictions for pp — > tiX^ on m^^, 
assuming ji — —200 GeV, M2 — 300 GeV and tan /3 = 30, which means that two chargino 
masses are about 182 GeV and 331 GeV, respectively, and mi^ increases from 342 GeV to 
683 GeV for m^^ in the range 100-600 GeV. One finds that the total cross sections for the 
^2X2 channel are always smallest and less than 3 fb, but the total cross sections for other 
channels are large and range between 10 fb and several hundred fb for most values of m^^. 
Especially for the iiXi channel, the total cross section can reach 1 pb for small values of m^^ 
(100 GeV < mi^ < 160 GeV), which is almost the same as ones of top quark and charged 
Higgs boson associated production at the LHC. However, when mi^ get larger and close to 
m^^) the total cross section for the t2Xi channel is the largest, as shown in Fig. 7. Moreover, 
Fig. 7 also shows that the NLO QCD corrections enhance the LO results significantly, which 
are in general a few ten percent. 

In Fig. 8 we show the dependence of the LO and NLO predictions for pp — > iiX^ on m ~- , 
assuming ji — —400 GeV, m^^ = 250 GeV and tan /3 = 30, which means that mi^ — AlA 
GeV and m-- increases from 417 GeV to 434 GeV for m-- in the range 100-300 GeV. One 
can sec that the total cross sections for all channels can exceed 10 fb. For small values of 
m^-, the iiXi channel has the largest total cross sections, and for large values of m^-, the 
tiX2 channel has the largest ones, which all can exceed 100 fb. Moreover, Fig. 8 also shows 
that the NLO QCD corrections enhance the LO results, and especially for the iiXi and t2Xi 
channels, the enhancement can exceed 50%. 

The associated production of ii and Xi is the most important since the total cross sections 
are the largest for m^^ < 400 GeV and m^- < 230 GeV. Thus we mainly discuss this channel 
below. 

Fig. 9 gives the dependence of the K factor (the ratio of the NLO total cross section 
over the LO one) on rrif^ for the iiXi production. Here we assume the same values of the 
MSSM parameters as in Fig.7. The K factor contains four contributions: (i) The curve 
(b) shows the ratio of the improved Born cross section, which is defined as convoluting the 
LO partonic cross section with the NLO PDFs, over the LO one. The corresponding K 
factor ranges from 1.18 to 1.39. (ii) The NLO cross section in the curve (c) contains only 
the contributions from the subprocess with two initial-state gluons, which are positive and 
the corresponding K factor varies from 0.30 to 0.06. (in) The NLO cross section in the 
curve (d) collects only the contributions from massless (anti)quark emission subprocesses as 
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shown in Eq. ()3.3U|) - ()3.32p . and the corresponding K factor keeps always about 0.10. (iv) 
The NLO cross section in the curve (e) includes only the virtual and the real gluon emission 
corrections, which are negative and the corresponding K factor ranges between —0.36 and 
—0.29. Finally, summing above four parts leads to the total K factor as shown in the curve 
(a), which indicates that the NLO QCD corrections enhance the LO total cross section and 
vary from 21% to 27%. 

In Fig. 10 we show the dependence of the total cross sections for the tiXi production on the 
renormalization/factorization scale, assuming /i = —200 GeV, M2 = 300 GeV, tan/3 = 30, 
mj^ = 250 GeV, and fir = f^f- One finds that the NLO QCD corrections reduce the 
dependence significantly. The cross sections vary by ±15% at LO but only by ±4% at NLO 
in the region 0.5 < fif/m.^^ < 2.0. Thus the reliability of the NLO predictions has been 
improved substantially. 

The cross sections of the iiXi production as a function of tan /3 are displayed for m^^ = 
200,300 and 400 GeV in Fig.ll, assuming /i = -200 GeV and M2 = 300 GeV. For the 
dotted curves, Arrib in Eq. fl4.3|l is replaced by Am'^ which only includes the corrections of 
the gluino piece and neglects the rest of the SUSY diagrams in Eq. fl4.4j) . From Fig.ll, one 
can see that the cross sections become large when tan j3 gets high or low, which is due to the 
fact that for the coupling b — ti — Xi low tan f3 the top quark contribution is enhanced 
while at high tan/5 the bottom quark contribution becomes large. Fig.ll also shows that 
the NLO QCD corrections enhance the LO total cross sections, and for m^^ = 200 and 
300 GeV, the enhancement is more significant for the medium values of tan /? than for the 
other values. Moreover, after comparing the solid curves with the dotted ones, we can see 
that the SUSY corrections to Anih except the gluino piece decrease the total cross sections 
significantly, especially for larger values of tan /?. 

Fig. 12 shows the dependence of the total cross sections for the tiXi production on the 
parameters /i and M2, assuming m^^ = 250 GeV and tan/5 = 30. Note that the phenomenol- 
ogy at CERN e~^e~ LEP and Tevatron [36] has excluded the parameter range |yu| ^ 180 GeV, 
where one chargino mass and one of the neutralino masses can become very small, and es- 
pecially the latter is zero for vanishing fi. So we focus on the range ^ 180 GeV. One can 
see that the cross sections increase with an increase of M2 for |/i| < 450 GeV, and decrease 
with an increase of \n\ except > 350 GeV for M2 = 300 GeV. The NLO QCD corrections 
can increase and decrease the LO total cross sections depending on the values of /i and M2. 
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In Fig. 13 we display the differential cross sections in the transverse momentum pt of Xi 
for the tiXi production in three cases, which are obtained after integration over all y in 
Eq. (|3.45|) . assuming /i = —200 GeV and m^-^ = 250 GeV. In all cases, we find that the NLO 
QCD corrections enhance the LO differential cross sections for low px, but decrease the LO 
results for high pt- 

In conclusion, we have calculated the NLO inclusive total cross sections for the associated 
production processes pp — > tiXk the MSSM at the LHC Our calculations show that the 
total cross sections for the tiXi production for the lighter top squark masses in the region 
100 GeV < m^^ < 160 GeV can reach 1 pb in the favorable parameter space allowed by 
the current precise experiments, and in other cases the total cross sections generally vary 
from 10 fb to several hundred fb except both m^-^ > 500 GeV and the ^2X2 production 
channel, which means that the LHC may produce abundant events of these processes, and 
it is very possible to discover these SUSY particles through the above processes in the 
future experiments, if the supersymmetry exists. Moreover, we find that the NLO QCD 
corrections in general enhance the LO total cross sections significantly, and vastly reduce 
the dependence of the total cross sections on the renormalization/factorization scale, which 
leads to increased confidence in predictions based on these results. 
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APPENDIX A 

In this appendix, we collect the explicit expressions of the nonzero form factors in Eq. ()3.5|1 
and Eq. ()3.(i|) . For simplicity, we introduce the following abbreviations for the Passarino- 
Veltman three-point integrals Cj(j) and four-point integrals -Dj(j), which are defined similar 
to Ref. [2^ except that we take internal masses squared as arguments: 

Cl(,) = aO)(0,0,s; 0,0,0), 

^iU) = ^^U) ("^1 . "^J- . 0, m| , 0) , 
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rte 

^4 ) ^ 

— 

rii — 
r~ik 



i,m|,0;m|.,0,m|.), 



m|-,0,t;m?,0,0), 
0,0,s;m?,m?,m|.), 
s, 0,0;m?,,m?,m|.), 

J. 2 n 2 2 2\ 

t,m^-.,0;m^,m(,m-), 

J. 2 n 2 2 2 \ 



^4) 



0') 



= D. 



i(3) 



0,mf.t;0,0,m.2 



2 n 2 2 2 \ 



m|,i,0,s,0,m?-;0,m|,0,0), 



ml , M, 0, s, 0, m ; 0, m? ,0,0), 



t, m^-,u, m^- , 0, 0; 0, m^^, 0, „.j^^. 



AJ.2 2n2222\ 

0, t, m-- , s, m^. , 0; m^, m^, , . ) , 



2 ri2n22 22\ 

s, m-- , M, 0, m^-^ , 0; m^, . , , , ) , 
t, ml-,u, ml, 0, 0; m?, m^, m?., m^^ 



Many above functions contain the soft infrared and coUinear singularities, but all Passarino- 
Veltman integrals can be reduced down to the scalar function Bq, Cq and Dq. Here we 
present the explicit expressions of Co and Do, which contain the singularities and are used 
in our calculations: 



^0 



C I 1 s 1 , s 27r\ 
— (— --In — + -ln — -^), 
s € mf 2 mf 3 



t 



Li2(-) + -ln' 



mf, ti 2 mf. 



2 



2 



^)-Li4)-lln^l 



3 1 S 1 "^f ~ 2 — ti TT^ t 

^[7-^ - - In — + - In — ' ^ - - In — - — - Li2(-) + Li2(^ 2" 

sti 2e^ e mf e mj e mf 2 ti mj — mi- 

H i^i t^i Xf^ 



-Li2(l + 



1 



mf — ml- t . 

_h ^^)+2Li2(-A)_iln2 *^ 

s ' ^ sti' 2 m~ 



9 22 
ml- „ m| — mt- 

2-^ + ln^ln^^^ 
I. "^1. 

Zi Zi Zi 
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+^ In" 3^ + ln "^\; "° (ln - In -^) - R(0, 1, -ti, t2, 0, -£2") 

,2 ^2 



mf sti s mf mf 



4- fTl~ — TTl 

-R(0, 1, St,, h, -1, --) - R(0, 1, St,, h, -1, 1 + ' 



1 S 

L>g = L>S ^ U, ti ^ lii, t2 ^ ^^2), 



t2 



^0 = -^[tt:^ - - In + - In -2-^ - 7 In —o- - Li2(-) + L12 



H H X/c 



?/ 1 —TTii- 1 1 „, _f m-_ 

-L12 — )--ln 2 — ^ + -ln — + -ln + 1, s, 1, f 

1*1 2 rrij 2 mi 2 mj mi. mi. mj 

222 

+ _f mf — mi- mi- 

+R( — , 1, U, 1, -1) + R(0, 1, mti, /e, ' 2 ""' ^ ^) 

-R(0, 1, -uiii, /c, ^) - R(0, 1, wi^i, /c, — -^)], 

where we define = r(l + e)(47r//r/?77.|,)^, and 

/b = ^(^1 + ui) + ml {ml - ml- ) , fc = m\m\- - tu, 

R{xi, X2, X3, X4, X5, xq) = Li2(^^^) - Li2(^^^) + In — ^ In ^ 

Za Za X4 

with A' = X3XQ — X5X4, x[ — xs + X1X4, and X2 — xs + X2X4. 

Moreover, there are the following relations between the form factors {i — 1...6): 

Thus we will only present the exphcit expressions of and f2°-i^^- 
For the diagram (a) in Fig. 2, we find 

n = ^4{(9 - 2e)So(s, 0, 0) - 40(1 - 2e)C^^ + 18sC« + (i?5i)'[36Co^o + 4Co^o 
-185o(s,m2,m?.)]}, 

/s = 2/" + ^4{«(7Cr - 2C« - 20(1 - e){Cl2 + C"^ + C|) + 16C|) + {R],f[2Cl{ml 

-ml) + 2i?o(s,m?,m|) - 25o(0, m|,, m|,) + 2s(Co^ + Cfi + Cfa + 2Cf + 9C| + Cf 

+9C7/2 + 9^^24)]}, 
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For the diagram (b) in Fig. 2, we find 

= ^{2l%R],[m,RUCi + Ci) - m^li,{Cl + Ci + Ci)] - m--[4(2Co^ + + 20^) 

fs = ^{4[5o(m|- , 0, m|) - 2So(s, 0, 0) - sa(2C', + C\ + + m| (C? - 2Cl + C^^) 

-m^.C^I + 2k%R%\mmAiCl + 4(i^o(m|-, m?, m|) + m^C^ + m? (C^ + C^))] 
-2(^2 + «2 + m^_)[4(2C^ + + 2C\) - 21^,1^,1^0^ + 2m^-j],R],[rn-A{Cl 

ft = Vl 

For the diagram (c) in Fig. 2, we find 

n = -^4[2(l - 0, 0) + 2{R),fB,{s, m?, m|)], 

/s = 2/3^. 

For the diagram [d) in Fig.2, we find 

ft = [165o(t, 0, m|) + I6S0K, 0, m|) + ti(18Co^ - 2(7^ + 210^) - 8m|(C? - OC^^)] 

-7^^Ni(^o(t, m?, ml) + Bo{m%ml ml)) - (i?*i(2m? + 2m? - t - m|) 

'-J 

-4m5mti?y (9Cf - m|, 0; m?, m?, m^))]} + (4 ^ i?*^, ^ i?*^), 

/r = -/5 ■ 
For the diagram (e) in Fig.2, we find 

f! = ^{4[^o(m^- , 0, m|) - h{C^, + CI + C|) - 2m|Co^ - 2m^- (C,^ + 

-24i?5JA;5,(i?o(m^-, m|,) + m^C^ + t^Q + m|-q) - /,>^-m,(C^ + Q)] 
+24i^im*(A;5,m^C^ - 4m-- Ci)}, 

/r = -/s ■ 

For the diagram (/) in Fig.2, we find 

O It Alt 

^ "3^^f [^^^ + °' ""^'^ " ^"^""^^^^ + T^^^^^' + ^^^""'^ + ^lBo){R!i,R], 

- mtm5(44i + 442)^0] {t, ml m])}, 
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For the box diagram (a) in Fig. 3, we find 

^Box(a) ^ i_4^._[5^(^_^ 0, m|) - Bo{t, 0, m|) + ^^(Ct + 4D«o - 2iiD«)], 

/s"""^^"^ = |-4{2So(m|-, 0, m|) - 2So(i, 0, m|) - i2[2Co" + 2Q + C^^ + 2(7^ - 2Q + C^^ 
+6L'«o + s{D1^ + SD^'g + 2D1 + 4D^3 + 4D^) - sa{D1^ + 2D?3 + 2D» + + 2D^) 
-h{2D^ + 3D?2 + 3D^ + 3D^3 + 2D^ + 3D^) + U2{2D^ + L'^i - SD^^ + ^^3 + S^*? 
-3L'«3 - 2i^2 + ^3)] - t2ml (2L>?i + L>«3 - 4L>« - L>3«3 - GL*") + ism^- (L'?3 + L>3"3 
-2^3)}, 

A^"^^"^ = -|^4{4i?o(m|- , 0, m|) - 4Bo(t, 0, m|) + t^P^o^ - 4Co^ + - s{4Dt, + 2D? 

+ SA(i^?3 + 2i^? - 2L'3« + L>?i) - - 2L>«) - ^^2(3L'? - 2L>« + D^,) 

-{2DI - 3Dt,){t^ + U2)] - hml{2Dt, - D^, - 4D? + 2^3^) + t^m^- {Dl, - 2DI)}, 

= 64m-- {Dl^ + Dl, + D»3 + D3«3 + Dl), 
/r^") = -64m,-p?, + i5?3). 

For the box diagram (h) in Fig.3, we find 

^Box(6) ^ _i4^._(2Co^ + d - AD',,), 

A"""^^'^ = ^4[2Co" + 2Co^ - Coin, , 0; 0, m|, 0) + + - AD',, + (t + m|-)L>^ 

-2uiD', - (Ml + 2m|)(Dj + Dl)], 
A"""^^'^ = ^4[2(^?o" + C^i + C^2 - D'oo) + + t2){2D', + D\, + + 3DJ + 20^ + D^) 

-s(L>^3 + Dl^) + + ^3) - 2m|(L>Ji + L>?2 + ^13 + ^1)], 

/r''"^^'^ = ^4[2^^o + Cl + m\- , 0; 0, m|, 0) + u,{D\, + 2^^^) + 2m|Dji], 

fir''' = -li\u^^-D\,. 

For the box diagram (c) in Fig.3, we find 

.Box(c) _ .Box(c) _ _2 t 

^Box(c) ^ _ 1 4|2So(m^ , 0, m|) - 2Bo{t, 0, m|) + i2[-2Co^ - 2Ct + 2D^,, + {s + 2t 
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-^3^3 + ^3^) + Dl^{t + mj, - 2u2) + Dl^{t + u + 7/2)]}, 

/r''"^'^ = 3^4{2i?o(m|- , 0, m|) - 2Bo{t, 0, m|) + i2[-2Co^ + sD^ + (t^ + u^) + i^J^ 
+£'^3) + (Co + Ci + C2)(mJ-, 0, ^; 0, m|, mf.) + iii(2L>^ + 3L>J + L'^i + 2L>J3 

+D^3 + 3^3^) + 2ml- + ^^^3 + D{ + D^^ + D^)]}, 

For the box diagram {d) in Fig. 3, we find 

+W2/^o' + "2/^2') + ^*4(C^o - 2/^00 - ^2/^2')]}, 
+mlDt + tDi)]}, 

+k%[mtm-,Rl,{Dt + Dt) + 4(C^' - 2Do'o - ^K^^ + ^2) + (^i + «2)(i^^i + ^^2 
+L>f3 + L>f + L>2'3) + - mliDf, + D^))]}, 

+m?D2^ - h{Di^ + ^2^3) - «2(i?^2 + ^23) + ^1^22)]}, 
^Box(d) ^ e^J^i^Sj^.^t^p^d ^ ^ 2L>^2 + L>^3 + 2L>^ + Di^ + L>2'3 + 2i^2') 

-mtRlM + ^Df, + Df, + Df + Di^ + D^^ + Di)] - k^m^-RUDt, 

+Dt, + Dt, + Dt + D',,)}, 
^Box(d) ^ 6i?5,{Z5Jm,4(L'f2 + ^22 + ^23) - m^RlM, + D^, + D^, + L*^')] 

For the box diagram (e) in Fig.3, we find 

9 - - 
pBox(e) _ ^7,b pb n*^ 
J3 — '^'^jk^jl^i2^mi 

^Box(e) ^ [(4^^ _ m~X,){Dl^ + i^i2 + + ^2) - ^§4(^0 + Dl 
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-rDl + Dl)] + 4[4(2i^oo + ^(^12 + ^22 + ^23 + ^2)) + «(^12 + ^11 
+1)^3 + + m^(4m5 - 4^t)(^o + ^1 + Dl + ^3^)]}, 

+Ci + C2) (m^- , 0, u] ml m| , m| ) - m?L>2^ - - m| D^s)]}, 

+m~,RUD^, + Dt + Dl + Dl)] - k%m--li,{Dl, + Df, + Df, + D^)}, 
^Box(e) ^ lRi^{i^^[m,Rl,{Dl, + Dl,) - m^RiiDl, + Dl, + Dl)] + k%m^-4,Dl,}. 

For the box diagram (/) in Fig. 3, we find 

^Box(/) ^ ^Ri^{ii^[(^rntRi, - m~,Ri)(C^ - 2DI, + s{dI + D( + D^ + D^)) - m,4(C^ 
+sDl + hDl)] + k%m^-R\,[Cl + 0, m^, m^, m?.) + - 2L»o^o + sDi 
-U{Di + D{ + Di)]}, 

^Box(/) ^ 1 i^5^|/5^^._ + d{ + D3O - m,4(L>{ + L'a^)] + ^[^(^^ - ^0 

-Co(0, mj- , ix; m^, ml m?,) + 2/^1(0 - ^I^o - sDi + + dI) - ml{D{ + Z)30) 

^Box(/) ^ [(^^4 _ m54)(L'{i + D{, + 2L>{3 + + ^4 + ^^3 + dI) 

-m-gRl.iDl + D( + D| + D()] - k%[m~,{m,Ri - m-^I^M + D{ + Di + D() 

+4(C^ + q + 0, m^-; m^, ^^'^ + - 2i?i5o + (s - + D{, 

+D(, + D() - ml{D(, + D(, + Di, + D^))]}, 
^Box(/) ^ 2^5^{/5^^_^_ [mgRl.iDl + D( + D{) + (m^i?*! - m,i?*2)(/^{i + 2^(3 

+D( + A(3 + Dl)] + k%[{m-,mtRl, - - h)){Dl + D( + D^) 

-4(-(5o(m|- , ml ml) - Bo{u, ml ml)) - + -C\- s{D{^ + D{, 

+D( - D() + m|(D(i + 21^(3 + D[ + ^{3 + Di))]}, 
^Box(/) ^ + D( + Di + D{) + (m^4 - m,4)(D(3 + + ^3^ 

+L>30] + k%m--R^,{D(, + D(, + D(, + D()}, 
^Box(/) ^ ^i?5^{/5J^,4(z^/3 -D( + Di,) - m-XM, + ^{3 + Di)] 

-k%m^-RUDi + D(, + D()}. 
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For the box diagram (g) in Fig. 3, we find 



fBox(g) _ 7 J 
J5 — o ik^'i- 



APPENDIX B 

In this Appendix, we collect the explicit expressions for the amplitudes squared of the 
radiations of a real gluon and a massless (anti-) quark. Since these expressions are only used 
to the hard non-collinear parts of real corrections in Eq. (|3.16p and Eq. (|3.35|) . they have no 
singularities and can be calculated in n = 4 dimensions. 

For the real gluon emission, we find 



|Mf;| =Xik[\M(')\ +|MW| + |M(^*)| ] + (s ^t', S4 ^ M7, S3A ^ U2), 
where Xik = + \Kk\^)^ and 

2 1 1 

= { [t2u'ml - SiUQ{t2 + U2) +M6(m6(s3A - ^^2) + (^5 + (-^SA + ^2) 

-m'(s5A + - S4t'ml-] - -^[s^At'uQ + Ug(u2 - S3a) - UqUjIs^a + h 

+U2) + (54^2 - ■5S5a)(u6 + M7) + t2Uj - m|(M2(2s + t') + At2t' - SsS^a) 

t' 

-ml- i2ue{t' + s)+ Ur{4:S + t'))] - -[At{4:t2 + U2) + (4^6 + n7)(2mj_ + ssa)] 

u' 1 

-[Att2 + M6(s5A + 2m|-)] - — — [IIS5AM6 - "^|(31t2 + 34^2) - m|-(9M6 



36s5t' '^fc 72s4 

+34U7)] - ;^[Ai(31t2 + I6M2) + (31^6 + 16m7)(s5a + am^.)]} 

+ ~2{^~79 [4'5(^t'53A + S4S5A + 2s4m|_) + u'At{u2 - 2^2) + tt'(s4 - 2u6)(s5A 



+2m^„) + S5m'(s5a - "^Ir)]} + ^"^1' 



'4 

72 1.1 



|M(*)| = ^{Y^[2t2(s3A - S4 + + ^2 - U7) + 4t2u'(s4 + S5A + ^6 " ^) + SsA^im" 

-T^K(S3A + S5A + + tte) - ml{s - S3A - t' + M2 + ^7)] + [t2{Am\- 

As^U ' 3054^7 ^k 

1 9 

+S4M7) + 2s5Atim|] - Y^[m|(i:2(16m| + Ts^) + 8ti(s3A + Sja)) - -S4S5Ati]} 

7[t2M7 + 4s3Am| + AujW? - S^a{s + S3A - t' + ^2 + M7) - M6(2S3A 



16S4tlM""^'' ■ ""-^^ ■ ^k 

26 



-2t2-U2)l 

~2 1 r 9 r / ., ~ 2 /j./ I I 2 



+(S5A + 2ml-){t' - s))] + -^[16ml{u2{t' + Uq) + t2{2t' + 2ue + U7 + u')) 
+16m|- (t'(2uQ + ur) + Uq{2uq + 2u-j + u')) — s4'2{Qs^a + 2uq — lu'j — 7m|) 

+ 7S4M2M6 + SS4(7S5A - 2m^-) - (9S5(s5At2 + 3t2«6 + 'W2M6 + ^2^7) + 1455^2^1 



+S4m^-(9t2 - 14ii6) - 2s'U6(s5A + 2m^-) - 2smf.(8s3A + 7s5A + 14m^-)] 



r2 

9 



^ -[2t'(Af + S5A)(2t2 + U2) + (S3A - S4 - 2^2 - - U2 + 2^6 + M7)(^2M7 
-W2ti6) + (S5^2 - 'Wm\-){2UQ + -Uy) - 3sS4i2 " 2sAt(2t2 + U2) — 2sS5A{t2 

9 

+U2 - Uq) + 4sm^-(s4 + 3^6) + (S5 - s)u2Uq] - ;^[sS5a(s3A " S4 + t') + (s 
+'U2)(S4^2 + SsaUg) + SAU2U6 - ^2^7) + 2i'm^- (s + 1^2 - U7)] [tlUeisQA 

SU-J 

+2ml-) - ul{u2 - 2ml-) + ^i^^m? + u^ml{2t2 + U2)]] + ^^{^[2«6(^2 

9 

+1*2) - S^aUqu' + t2ti'm|.] - — ^[(2^2 + U2){U(<,{t' - SzA + ti2 + U-j) + 2t'(s5A 

„ 16m? 
-2mf.)) - t2U'!{s3A + 4^2 + i + 3ii2 - 1*7)] + — T^[ui{t' + lie) + ^2(2^' + 2uq 

' 5411 
+1*7 + u')] - --^[2M6m|(2t2 + U2) + tl(s5AM6 + '"2"T'|.)]}- 

For the subprocess with two initial-state gluons, we find 



= -- |Mf^ I (s u', S4 ue, S3 ^ Si + uy + u', ssa ^ S4 + 1*7 + li' + At), 
o 

For the subprocess with the initial-state qq, we find 

jMfff = ^Xik{^[At{t'u2 + t2u') + {t'm + ^/6^i')(s5A + 2m|-)] - ^[SS3(A, + S3 

-S4 - S5a) + 53(^2 + - Uq){u-j -U2- u')] ^[s^At'u2 - 2sS3(At + S^^a) 

-54^2(2^2 + t' - ue) + t'u7(s3 - Ssa) - 531*7(2146 - ^2) + 531*6(1*2 + u') - 54^2(1*' 
-1*7) - 2ml-{t'{u7 - U2) - 2SS4) + u\t2 + 2t' - Uq){s^a + 2ml-)]}. 

^k ^k 

For the subprocess with the initial-state hq or hq, we find 

= VWhs ^ U\ 54 ^ U,, t2 ^ S3A, S3 ^t2,-^^). 

S3 tl 
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For the subprocess with the initial-state bb, we find 



-2 



1 ^ rS-SsA 



= l^Xik{^[s{Si + S5A - A, - S3a) + {t2 + t' - U6){U2 - U7 + u')] 
Z I s 

2t' 1 

H —[s5AU7 - U2m}, + (m2 + 2M7)m^-] -[^st'u2 - SS3AS5A 

+ S2,aUq{u2 - My) + 54^2^7 + (^SA " S'^a)'^^! + (^2 + ^){-SaU2 + Ssam' 

1 



Af(sS3A - t2y!) + 2(SS4 + (t' - M6)m')"^|-] + 1 "I 7[2sS3a(s3A - S4 

S5A + t') - t M2(2t2 - 2S4 + t') + (2t2M' + ^'(^2 - 2^7 + 4Mp)) + 2(^4^2 



6 

-■53Ai2)(-S ^U2- U-j) + S^A^{u-j - 2s) + m^i^t^VL - il U2)\ + ^-^[At(t'M2 + t2^i') 

3 

+ (t'M7 + M6W')(S5A + 2m^-)] [s5At'(u2 - U7) - 2sS3A(At + Ssa) 

-S4M2(2t2 + t') + SsAUrit' - 2Uq) + (S3A + S4)(m2M6 + ^2^7) " S4t2M' + S^au' {t2 

+2t' - Me) + s^au&u' + 2m2 (2ss4 + (^2 + 2t' - + t'{u2 - M7))]} 

+ (S ^ m', S4 ^ Me, ^2 ^ S3A, S3 ^ ti). 

For the subprocess with the initial-state 66, we find 



Wik\ (■« ^ -54 ^ M7, M2 ^ S3A, S3 ^ Ml). 
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FIG. 1: Leading order feynman diagrams for gb — > Ux^ ■ 
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FIG. 2: Virtual one-loop Feynman diagrams including self-energy and vertex corrections for gb 

UXk ■ 
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FIG. 4: Feynman diagrams for the real gluon emission contributions. 
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FIG. 6: Dependence of the total cross sections for the iix^ production at the LHC on the cutoff 
Ss, assuming = -200 GeV, M2 = 300 GeV, m^^ = 250 GeV, tan/3 = 30 and Sc = Ss/WO. 
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FIG. 7: Dependence of the total cross sections on rrif^ for the iiX^ productions at the LHC, 
assuming n = -200 GeV, M2 = 300 GeV and tan^ = 30. 
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FIG. 9. K = ctnlo/'^lo factor for the tiXi production at the LHC as a function of "i^^, assuming 
/i = —200 GeV, M2 = 300 GeV and tan/5 = 30. (Jnlo is the NLO total cross section for (a), the 
improved Born cross section for (6), the corrections of the subprocess with two initial-state gluons 
for (c), the corrections of the massless (anti)quark emission subprocess Ea. ()3.3U() - H3.32() for {d), 
and the virtual and real gluon emission corrections for (e). 
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FIG. 10: Dependence of the total cross sections for the tiXi production at the LHC on the 
renormahzation/factorization scale, assuming /j, = —200 GeV, M2 = 300 GeV, tan/3 = 30, rrif^ = 
250 GeV, fir = fif and rriav = (^nf^ + m^-)/2. 
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FIG. 13: Differential cross sections in the transverse momentum pr of Xi for the iiXi production 
at the LHC, assuming fi = -200 GeV and m^^ = 250 GcV. For (a), tan/3 = 30 and M2 = 600 
GeV; for (6), tan/3 = 30 and M2 = 300 GeV; for (c), tan/3 = 4 and M2 = 300 GeV. 
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